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Abstract 

After electron-positron annihilation and prior to photon decoupling the magnetic Reynolds num- 
ber is approximately twenty orders of magnitude larger than its kinetic counterpart which is, in 
turn, smaller than one. In this globally neutral system the large-scale inhomogeneities are provided 
by the spatial fluctuations of the scalar curvature. Owing to the analogy with the description of 
Markovian conducting fluids in the presence of acoustic fluctuations, the evolution equations of a 
putative magnetic field are averaged over the large-scale flow determined by curvature perturba- 
tions. General lessons are drawn on the typical diffusion scale of magnetic inhomogeneities. It is 
speculated that Reynolds numbers prior to electron-positron annihilation can be related to the en- 
tropy contained in the Hubble volume during the various stages of the evolution of the conducting 
plasma. 
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In a conducting plasma, such as the early Universe, the kinetic and magnetic Reynolds numbers 
are defined as [1, 2, 3] 

ry f rms Ly V^ms Lb „ -f^magn r-.\ 

-f^kin — ) -f^magn — ; ^magn — ^ i 

^th ^^magn -'I'kin 

where v^.-^^^ estimates the bulk velocity of the plasma while Vth and t'magn are the coefficients of thermal 
and magnetic diffusivity; and Lb are, respectively, the correlation scales of the velocity field and 
of the magnetic field. In Eq. (1) Prmagn denotes the so-called magnetic Prandtl number [1, 2, 3]. 

Prior to electron-positron annihilation (i.e. T > McV) the coefficient of thermal diffusivity can 
be estimated as z-'th ^ (a'^^T)^^ from the two-body scattering of relativistic species with significant 
momentum transfer. The conductivity of the plasma is cr ~ T/oem so that the magnetic diffusivity 
becomes t'magn = aem(47rr)^^. Assuming, for sake of simplicity, thermal and kinetic equilibrium of 
all relativistic species (which is not exactly the case for T ~ MeV) the kinetic Reynolds number turns 
out to be i?kin — O{10^^), the magnetic Reynolds number is -Rmagn — 47r/aguji?kin ~ 0(10^^) and 
-P^magn ~ 10^. The latter estimates have been obtained by assuming, in Eq. (1), Ly 2± Lb ~ 
(where is the Hubble radius at the corresponding epoch); when the evolution of the background 
geometry is decelerated (i.e. a{t) ~ t"^ with < e < 1) = the particle horizon coincides with the Hubble 
radius up to an immaterial numerical factor (i.e. e/ (1— e)) which shall be neglected throughout^. In the 
symmetric phase of the standard electroweak theory, picking up a temperature T > 100 GeV where all 
the species (including the Higgs boson and the top quark) are in thermal and kinetic equilibrium, z^th 
and t'magn can be computed [4] in terms of the hypercharge coupling constant since the non-screened 
vector modes at finite conductivity are associated with the hypercharge field. In the electroweak case 

iikin ~ C(10^^) and i?„agn ~ ©(lO^^). 

The hypothesis of primeval turbulence has been a recurrent theme since the first speculations on 
the origin of the light nuclear elements. The implications of turbulence for galaxy formation have been 
pointed out in the fifties by Von Weizsaker and Gamow [5] . They have been scrutinized in the sixties 
and early seventies by various authors [6] (see also [7, 8] and discussions therein). In the eighties it 
has been argued [9] that first-order phase transitions in the early Universe, if present, can provide a 
source of kinetic turbulence and, hopefully, the possibility of inverse cascades which could lead to an 
enhancement of the correlation scale of a putative large-scale magnetic field, as discussed in in [10, 11] 
(see also [12] and references therein). The limits i?kiii ^ 1 and i?magn S> 1 are customarily assumed 
in the scrutiny of hydromagnetic turbulence where both the magnetic flux and the magnetic helicity 



^In what follows we shall assume a conformally fiat Friedmann-Robertson- Walker background geometry "g^^ — 
o^(t)?7^i/ where ?7^i/ is the Minkowski metric, t the conformal time coordinate and a(T) the scale factor; the confor- 
mal time coordinate t is related to the cosmic time t as a{T)dT = dt. 



1 



are conserved since z^magn 1 i-e.^ 

j^B-dt = -fmagn ^ V X (V X • cZS, ^ d?xA ■ B = -2t'inagn d^x B{-V X B), (2) 

where V and E are a fiducial volume and a fiducial surface moving with the conducting fluid; B and 
A denote the comoving magnetic field and the comoving vector potential. In the ideal hydromagnetic 
limit (i.e. a oo, t'magn and i?magn oo) the flux is exactly conserved and the number of 
links and twists in the magnetic flux lines is also preserved by the time evolution. If i?kin ^ 1 and 
-Rmagn the System is still turbulent; however, since the total time derivative of the magnetic 

flux and of the magnetic helicity are both ©(fmagn) the terms at the right hand side of Eq. (2) cannot 
be neglected. Finally, if i?magn S> 1 and i?kin ^ 1 the fluid is not kinetically turbulent but but the 
magnetic flux is conserved. 

For T > Mev the Reynolds numbers can be viewed as a measure the entropy stored in a given 
Hubble volume. The entropy stored within the Hubble volume Vh ~ 47ri7~^(t)/3 is directly expressible 
in terms of the Reynolds numbers: 

Sh = UsH-^ = '^,^^RL, (3) 

where s is the entropy density of the plasma. The adiabatic expansion implies that the total entropy 
over a comoving volume is conserved. Conversely Sh, i-e. the entropy stored in the Hubble volume, 
increases as the plasma cools down since the temperature redshifts as a"^ (where a is the scale factor) 
but the Hubble volume typically increases faster than both during radiation (i.e. Vh ~ a^) and 
during matter (i.e. Vh ~ a^/^). Up to numerical factors Sh ~ 0(10^^) right before electron-positron 
annihilation and Sh 0(10^^) in the symmetric phase of the electroweak theory. The maximal 
entropy stored in the Hubble volume today'^ is obtained by integrating the entropy density of the 
Cosmic Microwave Background radiation (CMB in what follows) over the present value of the Hubble 
volume: 

S, = ^TTs.H,'' ^ 1.43 X 10^« = ^tt'T^, (4) 

where Ty = 2.725 K. In the standard lore the huge value of the entropy contained in the Hubble 
volume is the result of an appropriate theory of the initial conditions, since the adiabaticity condition 
can only be mildly violated after inflation and for a standard thermal history. 

According to Eqs. (3)-(4) it would be tempting to establish a causal connection between the Hubble 
entropy of the CMB and the largeness of the Reynolds numbers. While such a connection cannot be 
excluded for T > MeV, Eqs. (3) and (4), taken at face value, would imply that the kinetic and magnetic 



^In hydromagnetic turbulence it is customarily assumed that Prmagn — 1 while the flow is incompressible (i.e. the 
bulk velocity of the plasma is solenoidal) [1, 2]. This is not necessarily the case in the early Universe, as we shall see. 

''We do not consider here the possibility of a gravitational entropy associated with a Hubble screen approximately 
saturating the Hawking-Bekenstein bound and implying S'screen ~ H^^Mp ~ O(10^^°) » Sj. 
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Reynolds numbers must reach their maximum for T <^ MeV. The latter conclusion is incorrect insofar 
as the thermal diffusion coefficient sharply increases after electron-positron annihilation while the 
conductivity is only suppressed as ^Jrjim^ where T = aT is the comoving temperature and rrie 
is the electron mass. While prior to electron-positron annihilation i^kin S> 1 and -Rmagn S> 1, after 
e^-e~ annihilation .Rmagn is still very large, i2kin gets smaller than 1 and Prmagn sharply increases. 
Indeed, prior to last scattering, the thermal diffusivity is dominated by Thomson scattering and the 
concentration of the charge carriers is not of the order of the photon concentration (as for T > McV) 
but ten orders of magnitude smaller i.e. = i^rij where n-y is the comoving concentration of the 
photons and r/b = 0(10""^'^) is the ratio of the baryonic concentration to the photon concentration 
indirectly probed by big-bang nucleosynthesis and affecting the present abundances of light nuclear 
elements. The thermal diffusion coefficient is then given by: 

where the electron-photon mean free path A^e and the ratio between the baryonic matter density and 
the photon energy density i?b are defined as: 



hea{T)a^e' ^^^'^^ 4Va;^o/Uo7 z -h 1 V 0-02258 j ' 

the baryon matter density pb is the sum of the matter densities and of the ions and electrons; a^e is 
the electron-photon cross section; uj^q = /ig^bo and J7bo is the critical fraction of baryons. For the 
sake of simplicity we shall adopt, in the explicit estimates, the following fiducial set of parameters 

(Obo, J^cO,^^deO, ^0, ris, ere) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088), (7) 

which are determined from the WMAP 7yr data alone [13] in the light of the vanilla ACDM scenario. 
After electron-positron annihilation the conductivity given by binary collisions can be estimated as 



a(r) = ai \ -— , Ac(r) = = 1.105 x 10« , (8) 

aemV^ealnAc ^ ' 2e^\TrnJ V0.02258y ' ^' 

where ai = 9/(87r\/3) depends on the way multiple scattering is estimated and Ac is the argument of 
the Coulomb logarithm. 

Let us finally come to a more detailed estimate of the velocity field prior to last scattering. The bulk 
velocity is defined as the center of mass velocity of the positive and negative charge carriers present 
in the globally neutral plasma^. The customary assumptions of hydromagnetic turbulence imply a 
solenoidal bulk velocity field with Prmagn — [1, 2]. Conversely, prior to photon decoupling, the 



®If the charge carriers coincide with electrons and ions, denoting with mi and me the masses of the electrons and 
ions the bulk velocity of the plasma is defined as Vb = {meVe + miVi)/{me + mi). The center of mass velocity of the 
electron-ion system is often called baryon velocity. 
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bulk velocity of the plasma is not solenoidal (i.e. V • 7^ 0) and Prmagn 3> 1. This means that the 
plasma is compressible and the divergence of the bulk velocity is directly affected by the large-scale 
curvature fluctuations. Using Eq. (6) and assuming Vrms = 1 the kinetic and the magnetic Reynolds 
numbers can be estimated as i?kin — 0-03 and i?magn — 1.30 x 10^° for a typical last-scattering redshift 
z ~ 1090. The upper limit obtained on i?kin by assuming that Vrms coincides with the speed of light 
can be made more stringent since the large-scale flow, prior to last scattering, can be determined from 
the evolution equations of the baryon-photon system. 

For T < MeV and around last-scattering the differences between the velocities of the baryons 
and of the photons are quickly washed out because of the tight-coupling between ions, electrons and 
photons. Recalling that H = 5,- In a = aH, to lowest order in the tight-coupling approximation the 
truncated set of hydromagnetic equations reads [14] 

i2b + 1 -rtb + 1 Pha 4(i?b + 1) 



J ■ E 4 ^ 

d^Si, = 3drtp - V • iT^b H T, dr^i = "idrip - -V • V^b, 

Pha 3 



drB = Vx (iLb xB) + WnV^S + V X (^] - —^V X [(V xB)xB], (10) 

V / 47renc 

where S]^ and denote the density contrasts of the baryons and of the photons in the longitudinal 
gauge; (j) and denote, respectively, the (00) and (ii) fluctuations of the conformally flat background 
geometry adopted in the present discussion. The total Ohmic current J obeys an evolution equation 
which can be reduced to a consistency condition as in the Eckart approach to relativistic thermody- 
namics; such a relation is given by: 

J = a(E + v,^xB + ^-:^), (11) 

where = a^ne. The system of Eqs. (9)-(ll) is supplemented by the evolution equations of the 
curvature perturbations which have been written elsewhere in the longitudinal gauge and even in full 
gauge-invariant terms (see [14] and references therein). 

In previous studies [3, 14] the emphasis has been to see which are the effects of the large-scale 
magnetic fields on the scalar modes of the geometry. The hierarchy between i?kin and i?magn after 
electron-positron annihilation suggests the possibility of addressing also the complementary part of 
the problem, i.e. the effect of the large-scale flow on the evolution of the magnetic field. Neglecting 
the terms which are quadratic in the magnetic field as well as the thermoelectric term (i.e. the last two 
contributions in Eq. (11)) the large-scale flow can be determined by using a standard WKB analysis 
giving, in Fourier space, 

Uyb(fe,'r) = M7^(fc,T)sin[/crs(r)]e-'''/''d, ^ " |' 
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(12) 



where rs(r) and k^{T) denote, respectively, the sound horizon and the typical scale of diffusive damping; 
MTz{k,T) encodes the normalization inherited from the (adiabatic) curvature perturbations which are 
the only source of large-scale inhomogeneities in the vanilla ACDM scenario: 

-Mnik,T) = 3'/^cfiT) - 2l^7^(r)7^.(^), Tnir) = 1 - ^ f a\T')dT' . (13) 

As a result of Eqs. (9)-(10), the velocity field (12) is not solenoidal. This is situation differs from 
standard hydromagnetic turbulence and it is closer to the situation of acoustic turbulence (see e.g. 
[16], first paper). Following the standard conventions [13], the correlation function of curvature per- 
turbations in Fourier space is: 

{n.{mM) = —^^^\q + 'P)Vn{q), Vn{q) = AJ-] , (14) 

q \qp/ 

where, according to Eq. (7), = 0.963 and An = (2.43 ± 0.11) x 10"^; Qp = 0.002 Mpc^^ denotes 
the pivot scale at which the power spectrum of curvature perturbations is conventionally normalized. 
Given the intrinsic inhomogeneity of the large-scale velocity flow, it is natural to generalize the kinetic 
and the magnetic Reynolds numbers of Eq. (1) to Fourier space by keeping the dependence on 
the wavenumbers in the velocity, the dependence on the redshift in the diffusion coefficients and by 
choosing 2± = 1/A;phys(r) where /cphys = k{z + 1): 

Rkm{k,z) = J\fMnF{k,z), Prmagn{k, z) = J\fin3.gn G{k, z) 

•^■^'^ " 2 33/4m2i^o'^''"'"'^ 20C(3)ae„,lnAc(T)r,blTy' ' 

F{k,z) = Ho [1 - 2cl{z)]Tn{z, zeq) (^)"°"' '" ^'^^"^ , G{k,z) = j^^, (15) 

V^^p/ kc^l [z) + -"-J ' 

where C(3) = 1.202; introducing the function y{z,Zeq) = [-^1-1- (-^eq + + '^) — 1] the function 
T-jiiz, Zeq) can be expressed as: 

T(.. ^-(.^U^^ 2 [y{z, 2:eq) + 1] [3y\z, Zeq) + 15^(2:, Zeq) + 20] ) 

r7^(.,.eq) - (^ + 1)|1 - ^ [2 + y(.,.eq)]=^ /• ^''^ 

In Fig. 1 the kinetic and the magnetic Reynolds numbers are illustrated for the set of cosmological 
parameters of Eq. (7). The three curves in each plot correspond to three different wavenumbers. For 
even larger wavenumbers both quantities are oscillating as it can be argued from Eq. (15). 

The hierarchy between the kinetic and the magnetic Reynolds number defines naturally a pertur- 
bative scheme where the evolution equations of the magnetic field can be averaged over the large-scale 
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Figure 1: The kinetic Reynolds number (left plot) and the magnetic Reynolds number (right plot) 
for k = 0.0002 Mpc~^ (short dashed line), k = 0.002 Mpc~^ (long dashed line) and k = 0.004 Mpc"^ 
(full line). In both plots on the vertical axis the common logarithm of the corresponding quantity is 
illustrated. 



flow. Consider the magnetic diffusivity equation (10) and neglect all the terms which are of higher 
order in the magnetic field intensity. An iterative solution of Eq. (10) can then be constructed as 

B,ik,T) = Y.Bt'\k.r), Gk{y) = e-^'^'-'^^-y (17) 

n=0 

X (-m.ni(-ahn{qm+Pm)Va{q,Ti)B\^^{p,Ti), (18) 

where, for simplicity, t'magn is assumed to be constant in time. From Eq. (18) the first few terms of 
the recursion are Bf'\k,T), B^^\k,T) and B^'^\k,T). The term Bf'\k,T) = Qk{T)Bi(k) where Biik) 
parametrizes the initial stochastic magnetic field obeying flux conservation (see Eq. (2)): 

{Bi{k) Bj{k ')) = —Pij{k)PB{k) 5(3)(A: + k '), P^J{k) = ^ij " k h- (19) 

Since the curvature perturbations are distributed as in Eq. (14), the correlation function of the velocity 
for unequal times can be written as 

{v.,{q,T)v0y)) = ^U{q,\T-T'\)5^^\q + p), U{q,\T - t'\) = v{q) t,5{t - r'), (20) 

where, to avoid confusions with vector indices, the subscript 7b has been suppressed. The function 
v{q) appearing in Eq. (20) is 

v{q)=T,V{q), V{q)=Ml{q,n)^rn{q) sin^ [grs(n)]e-2'?V<?I, (21) 
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where r* denotes the last-scattering time and the correlation time Tc is the smallest time-scale when 
compared with other characteristic times arising in the problem. Because of the exponential sup- 
pression of the velocity correlation function for r > (where denotes the Silk time [15]), Tc 
approximately coincides with r^. The form of the correlator given in Eq. (20) is characteristic of 
Markovian conducting fluids [16, 17]. 

Denoting with Hi{k,T) the magnetic field averaged over the fluid flow, the terms containing an 
odd number of velocities will be zero while the correlators containing an even number of velocities do 
not vanish i.e. (5f"+')) = i/f = and (fif = i/f 0. So, for instance, (fif ^) = 
while {bP) = Hf^ is 

X / dnQkir - Ti) / dT2 Gpin - T2) {qm + Pm) (Qm' + P 
Jo JO 

X {Va'{q',r2)va{q,n)) Bv{p'). (22) 

After averaging the whole series of Eq. (17) term by term the obtained result can be resummed and 
written as: 

Hi{k,T) = {Bf\k,T)) + {Bf\k,T)) + {Bf{k,T)) + ... = e-'^'^-«"-S,(^). (23) 
where the magnetic diffusivity coefficient fmagn = 1/ (47r(7) has been renormalized as 

— Mli^{k,T^)rn{k) sin^ (A;/fc*)e-2*^'/'=i, (24) 

and k* = l/rs(r*). The averaging suggested here has been explored long ago in the related context 
of acoustic turbulence by Vainshtein and Zeldovich [16] (see also [17]). Prior to decoupling, however, 
both (v^) oc An ^ 1 and i?kin ^ 1- The resum indicated in Eq. (23) seems then to be more plausible 
in the present case than in the one of a kinetically turbulent plasma with strong inhomogeneities. If 
the Markovian approximation is relaxed the velocity correlator becomes 

{vi{q,Ti)vj{p,T2)) = qiqjv{q)r{q,Ti,T2)6'' '{q+P), v{q) = -^Vn{q)- (25) 

Assuming for sake of simplicity, Rhir*) <C 1 and Csb(T*) ~ l/-\/3 we have that 

riq,ri,T2) = ^jcos [qc^in - rs)] - cos [qc^in + r2)]|e-«'''*(-i+-2). (26) 

The results (24) and (26) lead to physically equivalent estimates and support the conclusion that the 
diffusivity wavenumber is smaller than expected from the usual arguments of the magnetic diffusivity 
scale at last scattering. In fact, ignoring the contribution of the bulk velocity in Eq. (10) the diffusivity 
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wavenumber can be roughly estimated as ka- ~ y/Ti^a. The standard estimate of k^ must be compared 
with the diffusivity scale arising from the effect of the large-scale flow: 

where cZa denotes the (comoving) angular diameter distance to last-scattering for the typical set of 
fiducial parameters of Eq. (7). The standard analysis based on the magnetic diffusivity scale would 
imply that all the modes k > k^ are diffused. The presence of large-scale flow implies that diffusion 
operates already for k > k^- This means that the correct diffusion scale to be considered prior to 
decoupling is not k^ ~ C>(10^^)Mpc~^ but, at most, k^ ~ 0(50) Mpc""^ which is closer to the Silk 
damping scale but qualitatively and quantitatively different. 

In summary, prior to electron-positron annihilation the large-scale (turbulent) flow can only be 
determined indirectly from the features of the various phase transitions. After electron-positron an- 
nihilation the flow can be inferred directly from the evolution of large-scale curvature perturbations 
imprinted in the CMB anisotropics. The hierarchy between the kinetic and magnetic Reynolds num- 
bers prior to last scattering suggests an effective description of the evolution of pre-decoupling magnetic 
which encompasses the conventional approach solely based on the conservation of the magnetic flux. 
On a more speculative ground the present considerations suggest that the Reynolds numbers are 
somehow related to the largeness of the Hubble entropy during the early stages of the evolution of the 
plasma. 
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